Classical double ionization of atoms in strong fields 
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Recent high resolution measurements of the momenta of two electrons emitted from Argon atoms 
in a strong laser field show a strong preference for the outgoing electrons to have similar momenta 
and to be ejected in the same direction along the field. We discuss the final stages of this process 
in a classical model in the spirit of Wanniers approach to double ionization after electron impact. 
Starting from an intermediate state of high but negative energy the field opens up a channel through 
which both electrons can escape. Near the threshold for this process Coulomb repulsion favors an 
escape with the electron momenta and positions symmetric with respect to the electric field axis. 
Classical trajectory simulations within this symmetry subspace account for most features of the 
observed momentum distribution. 
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Double ionization of neutral atoms (in particular He, 
Ne and Ar) in the presence of strong fields has attracted 
considerable attention after it was noted that the ob- 
served yield is much higher than can be expected on the 
basis of an independent electron model It was pro- 

posed early on that correlations between the electrons 
should be responsible for this enhancement. A striking 
demonstration of such correlations is provided by recent 
high resolution experiments on the distribution of the 
ion-recoil momenta and the electron momenta ||] 

which clearly show a preference for the symmetric escape 
of both electrons to the same side of the nucleus. This is 
very different from the process in the absence of a field 
(Wannier- ionization, where the dominant ioniza- 

tion path has both electrons escape symmetrically placed 
on opposite sides of the nucleus. In particular, while the 
Wannier mode minimizes the electron repulsion, a sym- 
metric escape to the same side has considerable residual 
energy in the repulsion between electrons. So why is this 
channel favoured? 

It is by now generally accepted that the ionization 
takes place in two steps: a high excitation of one elec- 
tron followed by the double ionization event |p|-p^ . The 
highly excited electron is accelerated by the field and 
driven back to the core where it collides with the other 
electron, and transfers enough energy and momentum so 
that both electrons can escape from the nucleus. It is 
this rescattering of the electron that enhances the en- 
ergy transfer and thus the ionization rate. But it also 
implies that during the ionization process the interaction 
between the electrons cannot be ignored. 

The amount of energy the excited electron can gain 
during half a period of the electric field has been esti- 
mated in 1^: for a field intensity of 2.9 • IQ^'^W/cm? the 
estimate gives a maximal energy transfer of about QQeV . 
This is insufRcient to ionize the second electron directly, 
but as the experiments show, double ionization is possible 
nevertheless. This implies that there must be a mecha- 
nism by which the electrons can draw additional energy 



from the field so that asymptotically, once the pulse has 
ceased, the total energy of the system is positive. 

The related problem of double ionization in collisions 
or single photon excitation was studied in a classic paper 
by Wannier . He assumed that during the first stages 
a high energy complex of electrons close to the nucleus is 
formed from which then the ionized electrons escape. If 
the energy is close to threshold, they cannot afford to put 
energy into the mutual repulsion and the escape is with 
both electrons on opposite sides of the nucleus. More- 
over, he argued that their distance has to be the same, 
for any difference in position and energy would be am- 
plified, pushing the configuration towards single electron 
ionization. 

Essential elements of this discussion also apply to the 
present situation. At the end of the first step the rescat- 
tering of one electron to the core produces a highly ex- 
cited two electron complex. The estimates show [|j that 
the total energy in the system is insufhent to doubleion- 
ize immediately. However, if the electric field during the 
collision is non-zero, a saddle opens through which the 
electrons can escape. Because of their mutual repulsion 
the electron that gets to the saddle first has an advan- 
tage: it can cross the saddle while pushing the other back 
to the nucleus. Since energy is scarce near threshold this 
results in either single ionization or in another rescatter- 
ing event, but most likely not in double ionization. This 
suggests that the dominant path leading to double ioniza- 
tion has both electrons cross the saddle side by side. This 
singular process then acquires a finite probability when 
trajectories asymptotic to this configuration are taken 
into account. 

Therefore, we propose that near the threshold for dou- 
ble ionization the only path leading to double ionization 
has both electrons escape symmetric with respect to the 
field axis. With the field pointing along the x-axis and 
the electrons confined to the plane z — their coordi- 
nates are [x, y, 0) and [x, —y, 0) in position and {px,Py, 0) 
and {px, —Py, 0) in momenta. In this geometry the obser- 
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vation that both electrons escape to the same side is built 
from the outset. As an aside we note that this symmetry 
plane also contains the Wannier orbit, for a; = 0, as the 
symmetric escape perpendicular to the field axis. 

The classical Hamilton function for this geometry then 
is (in atomic units, with infinitely heavy nucleus and in 
dipole approximation) 



H{px,Py,x,y,t) =pI+pI + V{x,y,t) 



(1) 



with potential energy 



V{x, y, t) = ^ +^ + 2Fxf{t) cos{wt + . 

^x^ + 2/2 2y 



and the pulse shape 



(2) 



(3) 



where the duration of the pulse is taken to be four field 
cycles, Td = Stt/w. The frequency is uj/2 = 0.057 a. u., 
and corresponds to the experimental situation 1^-^ . The 
rescattering of the electrons leads to a highly excited 
complex of total energy E which every now and then 
is close to the symmetric configuration described by the 
Hamiltonian (1). Any configuration on this energy shell 
(for some fixed time t) as well as any phase (/) of the field 
is equally likely, and the experimental observations are 
averages over initial conditions and phases. 

As mentioned, for the weakest fields where double ion- 
ization is observed the rescattered electron does not bring 
in enough energy for double ionization. However, if the 
collision happens near a time where the field is strong, 
the electric field distorts the potential and opens a path 
for escape in down-field direction. This process can be 
discussed adiabatically for fixed external field since the 
motion of the electrons near the nucleus is much faster 
than the change in the field. The ionization can thus be 
discussed in the potential (2) with fixed field. Equipo- 
tential lines for the potential (2) at a maximum of the 
field for F = 0.137 a. u., corresponding to an intensity of 
6.6 • are shown in Fig. 1. The saddle is lo- 

cated along the line x = rs cos and y = rs sin 9 with 
9 = tt/6 or 57r/6 and at a distance r| = V^/Fmax where 
Fmax — maxt|F/(t) cos{u!t + (j>)\. The energy of the 
saddle is 



Vs = -6jF„,ax/V3. 



(4) 



For the extremal fields in a pulse of the above mentioned 
intensity this gives Vs — — 1.69a.u., so that within the 
adiabatic picture the saddle can be reached if the return- 
ing electron brings in at least 1.22 a. u. in energy. During 
a field cycle the saddle moves in from infinity along the 
line at 9 — 7r/6, moves out again to infinity after half a 
period and then moves in and out again along the line 
9 = 5tt/6 during the second half of the cycle. Ionization 
is most likely when the saddle is closest to te nucleus. 



A typical trajectory within the symmetric configura- 
tion is shown in Fig. g. During the ramping of the field 
the electronic motion is little influenced by the electric 
field, but during the third half cycle of the field the sad- 
dle is close enough to the electron orbits and ionization 
takes place. Once on the other side of the saddle, the 
electrons rapidly gain en ergy. T he saddle thus provides 
a kind of transition state 18| for the double ionization 
process: once the electrons cross it, they are accelerated 
by the field and pulled further away, making a return 
rather unlikely. Moreover, they can aquire the missing 
energy so that both electrons can escape even when the 
field vanishes. The field thus plays a double role in de- 
termining a threshold for this process: during the first 
stages of the rescattering process it provides the energy 
for the collision complex and during the final stages it 
opens the path for double escape. 

In the experiments it is not possible to monitor all 
details and intermediate stages of the process. Most in- 
formation is extracted from the distribution of final mo- 
menta p„ of the nucleus and pi of the electrons, where 
Pi + P2 ~ ~Pn Because of the symmetry assump- 

tion in the model the components perpendicular to the 
field vanish. The distribution of the parallel components 
can be calculated by averaging over all initial conditions 
of prescribed energy and all phases of the field. 

Classical scaling of the Hamiltonian (1) implies that 
the results do not depend on the initial energy E and field 
strength independently, but on the combination F/E^ 
only. The field strength F is set by the intensity of the 
laser. The initial energy of the two-electron complex is 
determined by the field dependent efficieny of the sin- 
gle electron excitation step and thus not directly acces- 
sible, although it can be estimated as in We there- 
fore fix F and vary initial energy. The results for fixed 
field strength F = 0.137 a. m. and an initial energy of 
E = — 0.58 a.u. are compared to the experimental dis- 
tribution in Fig. ^. The final distribution of momenta 
clearly shows the double hump structure indicating a 
preference for ionization parallel and antiparallel to the 
field. The maxima of the experimental distribution are at 
about p|| = ±1.5 a.u., whereas the numerical ones within 
the symmetric subspace lie at about ±1 a.u.. We take this 
close agreement as strong indication that double ioniza- 
tion can only occur in the neighborhood of the symmetric 
process discussed here. 

For lower energy, E = —1.3 a.u. the minimum at = 
almost vanishes and only a single maximum shows up 
(Fig. IJa). This corresponds to the experimental situation 
of a weaker pulse which evidently transfers less energy 
to the rescattered electron Within the symmetric 
subspace, there are two reasons for this change: at this 
lower initial energy the electrons have little kinetic energy 
when crossing the barrier, so that the splitting should be 
expected to be small. Secondly, the electrons cross the 
barrier typically when the field is strongest, so that after 
the ionization there is still considerable smearing of the 
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distribution due to the interaction with the remainder 
of the pulse. If the distribution is monitored immedi- 
ately after the crossing of the barrier the two preferred 
momenta parallel and anti-parallel to the field stand out 
clearly, as demonstrated in Fig. 

These numerical results in the reduced symmetry sub- 
space are in surprising agreement with the experimen- 
tal data. They show that the configuration with ro- 
tational symmetry around the field axis dominates the 
cross section. The interaction with the field shortly af- 
ter ionization is responsible for most of the smearing 
of the final distribution, additional contributions come 
from trajectories that are not symmetric. Actually, non- 
symmetric configurations are needed in order to obtain a 
non-vanishing cross section to begin with: the symmet- 
ric subspace is a set of measure zero in the phase space 
of the six degree of freedom system and acquires a finite 
overlap with initial conditions only due to non-symmetric 
initial conditions that are asymptotic to the symmetric 
subspace. The cross section for double ionization shows 
that this overlap is small ■ 

The picture proposed here for multiphoton double ion- 
ization is very similar to that of Wannier for double ion- 
ization through electron impact. The main difference is 
that now the symmetry between the outgoing electrons is 
a rotation around the field axis whereas it is a point sym- 
metry with the nucleus in the center in Wanniers case. 

Finally, we would like to mention that besides the dif- 
ferent symmetry there is another difference to the Wan- 
nier double ionization without field: In the zero field ion- 
ization there is only a single trajectory so that there in 
leading order semiclassical approximation no quantum 
interferences can be expected The present prob- 

lem falls into the category of quantum chaotic scattering 
|pO| , pT| where the classical ionization dynamics is chaotic 
and quantum interference effects between different paths 
cannot be ruled out. It will be interesting to pinpoint 
quantum interference effects in this system. 
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FIG. 1. Adiabatic potential V{x,y,t) for fixed time t in 
the symmetric subspace. The saddle moves along the dashed 
line when the electric field points in the positive a;-direction 
and along a second obtained by refiection on a; = during the 
other half of the field cycle. 



20.0 

15.0 

^ 10.0 

5.0 

0.0 
7 2.5 
3. 00 
LU -2.5 
-5.0 



(a) 



(b) 




0.0 



0.5 



1.5 



2.0 



1.0 

t/(27l/C0) 

FIG. 2. A typical trajectory in the symmetric subspace 
with E — —1.3 a. It. (a) distance of the electrons to the nu- 
cleus. The dashed line indicates the distance of the saddle. 
Note that before the double ionization occurs the effect of the 
field on the electrons is minimal, supporting the adiabatic 
assumption, (b) energy of the electrons. Note that the ini- 
tial state has negative total energy and cannot lead to double 
ionization. The energy increases once the electrons have es- 
caped from the nucleus far enough so that acceleration by the 
electric field dominates. 
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FIG. 3. Final distribution of ion momenta parallel to the 
field for F = 0.136 a. m.: (a) experimental distribution from 
(b) distribution from symmetric subspace with initial en- 
ergy E = —0.58 a. It.. The classical distribution is based on 
an ensemble of 2 ■ 10^ trajectories. 
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FIG. 4. Final distribution of momenta parallel to the field 
for F = 0.136 a.u. and initial energy E — — 1.3a.w.. (a) at the 
end of the pulse; (b) at the first zero of the field after cross- 
ing the saddle. This figure demonstrates the smearing of the 
distribution due to the final post ionization interaction with 
the electric field. It is based on 50.000 classical trajectories. 
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